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ABSTRACT 

Very recently, Verlinde considered a theory in which space is emergent through a holographic 

scenario, and proposed that gravity can be explained as an entropic force caused by changes in the 

^^ ' information associated with the positions of material bodies. Then, motivated by the Debye model 

Cn ^ in thermodynamics which is very successful in very low temperatures, Gao modified the entropic 

{^JQ. force scenario. The modified entropic force (MEF) model is in fact a modified gravity model, and the 

i~{ ' universe can be accelerated without dark energy. In the present work, we consider the cosmological 

constraints on the MEF model, and successfully constrain the model parameters to a narrow range. 

We also discuss many other issues of the MEF model. In particular, we clearly reveal the implicit 

VO ' root to accelerate the universe in the MEF model. 
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I. INTRODUCTION 

Very recently, Verlinde [1| considered a theory in which space is emergent through a holographic 
scenario, and proposed that gravity can be explained as an entropic force caused by changes in the 
information associated with the positions of material bodies. In this scenario, Verlinde has successfully 
derived the Newton's law of gravitation, the Einstein equations, and the law of inertia, from the entropic 
point of view. In fact, the entropic force scenario is similar to the old idea of Jacobson [2|, but also 
beyond it in some sense. Similar entropic insight into gravity has also been made by Padmanabhan [3[ 
independently and simultaneously. 

Here we briefly mention some key points of the entropic force scenario following [l|. Motivated by 
Bekenstein's argument [j , Verlinde postulated that the change in entropy near the holographic screen is 
linear in the displacement Ax, namely, 

Tn r 
AS^2TTkB^Ax, (1) 

n 

where m is the mass of test particle, whereas ks, c and h are Boltzmann constant, speed of light and the 
reduced Planck constant, respectively. The effective entropic force acting on the test particle due to the 
change in entropy obeys the first law of thermodynamics 

FAx = TAS , (2) 

where T is the temperature. If one takes the Unruh temperature T experienced by an observer in an 
accelerated frame whose acceleration is a, i.e., 

, „ I ha 

kBT=- , 3 

ZTT C 

to be the temperature associated with the bits on the holographic screen, from Eqs. ([1]) — 1^, it is easy 
to recover the second law of Newton 

F = ma . (4) 

Considering a sphere as the holographic screen, Verlinde assumed that the number of used bits on the 
holographic screen N is proportional to the area A = Airr'^, i.e.. 

According to the equipartition law of energy, the total energy inside the screen is 

E = ^NkBT. (6) 

Of course, one can identifies E with the mass M inside the screen through 

E = Mc^. (7) 

From Eqs. ((U, (H]), and ([S]) — (O, one can recover the Newton's law of gravitation 

F = G^. (8) 

where G can be identified with the Newton constant now. From Eqs. ([3]), (Hj) and (|8l), it is easy to find 
the gravitational acceleration 



GM 



(9) 



and the temperature 

T^-^f. (10) 

ksc 2tt ^ ' 

As shown in [1], a relativistic generahzation of the presented arguments directly leads to the Einstein 
equations. We strongly refer to the original paper [Ij] for great details. 

Soon after Verlinde's proposal of entropic force, many relevant works appeared. For examples, Cai, 
Cao and Ohta [5j , Shu and Gong [y] derived the Friedmann equations from entropic force simultaneously. 
Smolin Q derived the Newtonian gravity in loop quantum gravity. Li and Wang jSj showed that the 
holographic dark energy can arise in the entropic force scenario. Easson, Frampton and Smoot 
considered the entropic accelerating universe and the entropic inflation. Tian and Wu ^lOj . Myung 
discussed the thermodynamics of black holes in the entropic force scenario. Vancea and Santos 
considered the uncertainty principle from the point of view of entropic force. Zhang, Gong and Zhu [Ic 
Sheykhi il4il derived the modified Friedmann equation from the corrected entropy. Also, Modesto and 
Randono [13 discussed the corrections to Newton's law from the corrected entropy. Cai, Liu and Li |16| 
considered a unified model of infiation and late-time acceleration in the entropic force scenario. For other 
relevant works to entropic force, we refer to e.g. [17- 19, 36] and references therein. 

The works mentioned above are in fact closely following Verlinde's proposal of entropic force [l| . To 
be honest, here we should also mention the other works which are strongly criticizing the entropic force 
scenario. For instance, the author of [Sg] argued that there are some possible flaws in Verlinde's idea. 
In [33|, Culetu argued that the relativistic Unruh temperature cannot be associated with the bits on the 
screen in the form considered by Verlinde. In [40|, Hossenfelder argued that some additional assumptions 
made by Verlinde are unnecessary and there are some gaps in Verlinde's arguments. In |41| . Myung 
found that entropic force does not always imply the Newtonian force law, and the connection between 
Newtonian cosmology and entropic force cannot be confirmed. In j42i] . Li and Pang found that infiation 
is inconsistent with the entropic force scenario. In [43|, Lee argued that there are some inconsistencies in 
Verlinde's arguments from a classical point of view. 

So far, we have briefly surveyed the current status of the works relevant to the entropic force scenario. 
It is fair to say that the entropic force scenario is still in controversy. On the other hand, there is no 
breakthrough on entropic force after Verlinde's proposal [l|. A deep insight is needed to understand the 
nature of gravity. In addition, further discussions on the entropic force scenario are also desirable. Only 
when more and more results on entropic force are available, one can say something conclusively at that 
time. To this end, we would like to contribute our effort and try to learn more about the entropic force 
scenario. In this work, we will consider a modified entropic force scenario proposed by Gao [20|, which 
has some interesting features. And then, we will constrain the modified entropic force scenario with the 
latest observational data. 

In [201 1 , Gao noted that statistical thermodynamics reveals the equipartition law of energy does not 
hold in the very low temperatures. Instead, as is well known, the Debye model [2l|, [221 i^ very successful 
in explaining the experimental results when the temperatures are very low. Since the equipartition law of 
energy plays an important role in the derivation of entropic force, the entropic force should be modified for 
the very weak gravitational fields which correspond to very low temperatures. Especially, the large-scale 
universe is in such an extreme weak gravitational field, and hence the modified entropic force (MEF) 
makes sense in cosmology. 

Following [20] , we briefiy mention the key points of MEF model. Similar to the Debye model [2l|, [22] 
in thermodynamics, one can modify the equipartition law of energy in Eq. ^ to 

E = ^NksTDix) , (11) 

where D{x) is the Debye function which is defined by 

_y 

and X is related to the temperature T as 



^(^)".^^'^^' (^2) 



^^^, (13) 



in which Td is the Debye temperature. By definition, x is positive. With the modified equipartition law 
of energy, namely Eq. PT|) . similar to the original entropic force, one can easily to obtain 20] 



in which (nb. Eq. (fTO]) ') 



— ^ (14) 

r2 D{x) ' ^ ' 



^ = f- c^) 



where gn = (27r ks c/fi) Tn is the Debye acceleration. Actually, Eq. ([T^ corresponds to the modified 
Newtonian law of gravity. In the limit of strong gravitational field, g ^ go and hence a: <?C 1, from 
Eq. (fT^ it is easy to find that D{x) — ?> 1 and the Newtonian gravity is recovered. On the other hand, 
in the limit of weak gravitational field, g ^ go and hence x ^ 1, one can see that Djx) — >■ 7r''/(5x"^) 
and then g oc 1/^/r, which significantly deviates from the familiar inverse square law [20[. However, as 
argued in [20|, one need not to worry about the possibility of MEF against the experimental results of 
the inverse square law. Since these experiments testing the inverse square law were done on the Earth or 
in the solar system, which are actually in the strong gravitational fields, we have a; <C 1 and D{x) — > 1, 
therefore the deviation from the inverse square law are extremely tiny. The significant deviation from 
the inverse square law can only occur in the very large scale in the universe where the gravitational fields 
are very weak, and hence it can escape the detection of these experiments testing the inverse square law. 
Of course, this argument relies on a small g^. We will justify it later in the present work. 

Next, we turn to the cosmological issues in the MEF Model. For convenience, we set the units ks — 
c — h — 1 hereafter. Using the derivation method in p, [6|, one can find that the modified Raychaudhuri 
equation is given by [20| 



iTTG{p+p) = -[H-^ 



-2D{x) + -^ 



(16) 



where p and p are the total energy density and total pressure of cosmic fluids, respectively; K is the 
spatial curvature of the universe; H = a/a is the Hubble parameter; a = (1 + z)~^ is the scale factor 
(we have set oq = 1); 2 is the redshift; a dot denotes the derivatives with respect to cosmic time f; the 
subscript "0" indicates the present value of the corresponding quantity. Taking into account the Hawking 
temperature T for the universe p.2>\ 



from Eqs. (flBl) and (|T5]) . it is easy to see that 



T=^, (17) 



-=f^, (18) 

where iJ/j = f/^i. On the other hand, the energy conservation equation still holds in the MEF model, 
namely 

p + 3i7(p + p)=0. (19) 

From Eqs. p6p and (|19p . one can derive the corresponding Friedmann equation. It is anticipated that 
Friedmann equation is also modified, H'^ 7^ (87rGp)/3, due to the correction term —2D{x) + ?>x/{e^ — 1) 
in Eq. (fT6|. The MEF model is in fact a modified gravity model. Gao showed that the MEF model can 
describe the accelerating universe without dark energy. We refer to the original paper [20| for details. 

Since Gao has not considered the constraints on the MEF model in [20|, we will try to obtain the 
cosmological constraints with the latest observational data in the next section. Further, we will discuss 
some relevant issues of the MEF model in Sec. IIIII Finally, we give the brief conclusion and some 
meaningful remarks in Sec. IIVI 




FIG. 1: The 68.3% and 95.4% confidence level contours in the Q,mO — C parameter space. The best-fit parameters 
are also indicated by a black solid point. This result is obtained by using the data of 557 Union2 SNIa alone. 



II. COSMOLOGICAL CONSTRAINTS ON THE MEF MODEL 



In this section, we consider the cosmological constraints on the MEF model. To this end, we firstly 
rewrite the equations to suitable forms. Notice that we consider a spatially flat universe (namely K = 0) 
throughout this work. From Eqs. ([TO)) and p^ . we have 



8TTGdp = 3 



-2D{x) 



3x 



dH\ 



(20) 



As in [23, we consider the universe contains only pressureless matter. So, we have p = Pm = Pmoo, ^ 
Pmo (1 + z)^- Dividing 3Hq in both sides of Eq. (pO)) . we obtain 



where 



f^. 



rimo da 



-zhT 



^2D{x) 



3x 



E 



H 



c 

X 



dE^ 



C- 



Ho 



(21) 



(22) 



Note that Qmo =/= Pmo/pa ~ 1; because Friedmann equation has been modified in the MEF model, i.e., 
H^ 7^ 8ttGp/3. By definition, ( is positive. Finally, we get the differential equation for E{z), namely 



-2D 



3C/E 



;C/S-1 



• 2£; — = 3a„o (1 + z? 
dz 



(23) 



In principle, one can numerically find E{z) from this exact differential equation, and then fit it to the 
observational data to get the constraints on the MEF model. However, we find that it consumes a large 
amount of time beyond normal patience when we scan the grid points in the parameter space, mainly 
due to the hardness of numerically solving the exact differential equation (P5)) in which D{C,/E) is an 
integral whose upper limit is C/E itself. Therefore, it is advisable to find a reliable approximation of the 



exact differential equation ([23)) . Note that 

- 2Dix) + -^ ^l-K + O ix% (24) 

e^ — 1 4 

for any small quantity x. On the other hand, in [20| one might find the hint of a small ( (notice that 
it has been chosen to be 10~^ in [20| for example). Together with the well-known fact that usually 
E{z) increases rapidly when z increases, (/E is a small quantity in Eq. ([231) ■ Therefore, we find that an 
approximation of the exact differential equation (j23p is given by 

1- --^] ■2EdE ^flmoda^^. (25) 

4 E J 



Integrating Eq. (P5|) . we have 



3 

E'^ - -CE:^n„Ma^^ + const., (26) 



where const, is the integral constant, which can be determined by requiring E{z = 0) = 1. Finally, we 
find that 

E^^^CE = 0„,o {l + zf+(l-^C- a„o) , (27) 

which is a quadratic equation of E in fact. Noting that E is positive, we solve Eq. 1^7^ to get 



i^w = |c + i{^c' + 4 



f7„o {i + zf+(i--c-n„ 



1/2 

(28) 



Obviously, when C ^ 1, we see that the MEF model reduces to the familiar ACDM model in which 
E{z) = [Q,no (1 + z)^ + (1 - ^mo)] ■ Therefore, it is not surprising that in [20] Gao found the MEF 
model with C = 10"^ is degenerate to ACDM model. In fact, this observation is indeed the key point to 
understand the reason for accelerating the universe without dark energy in the MEF model. 

In the following, we consider the cosmological constraints on the MEF model from observational data. 
At first, we use the observational data of Type la supcrnovae (SNIa) alone. Recently, the Supernova 
Cosmology Project (SCP) collaboration released their Union2 compilation which consists of 557 SNIa [24| . 
The Union2 compilation is the largest published and s pec troscopically confirmed SNIa sample to date. 
The data points of the 557 Union2 SNIa compiled in [24'| are given in terms of the distance modulus 
A*obs(-Zi)- On the other hand, the theoretical distance modulus is defined as 

Mtft('Zi) = 51ogioL'L(zj) + ^0, (29) 

where /zq = 42.38 — 51og]^Q h and h is the Hubble constant Hq in imits of 100 km/s/Mpc, whereas 

^.W = (l + ^)f^, (30) 

in which p denotes the model parameters. The x^ from 557 Union2 SNIa is given by 

i ' 

where a is the corresponding la error. The parameter /iq is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over /iq. However, there is an alternative 
way. Following [25l |26| , the minimization with respect to /iq can be made by expanding the x^ of Eq. (j3ip 
with respect to /iq as 

xl{Y>) = A-2^i,B + ^ilC, (32) 
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FIG. 2: The same as in Fig. \T\ except that this result is obtained by using the combined observational data of 
557 Union2 SNIa and the distance parameter A from LSS. 



where 



Mp) - E 



[^J'obs{zi) - iJ,th{zi] Ho = 0, p)]^ 



B(p) = ^ ^^2^iM. 



fJ.thizi;Ho = 0,p) 



,(^0 



C = 






M) 



Eq. (15^ has a minimum for /iq = B/C at 



X^(P)-^(P)- 



C 



(33) 



Since xfi min — Xu min obviously, we can instead minimize xfi which is independent of /xq. 



model parameters are determined by minimizing the total y^ . When SNIa is used alone, we have x^ 



A/i 



which is given in Eq. (gSl)- As in [27, 2^ 
1.0, 2.3 and 3.53 for Up — 1, 2 and 3. 



J, the 68.3% confidence level is determined by Ax^ 
respectively, where Up 



X 



Amzn 



is the number of free model parameters. 
2 - -2 < 4_o, 6.17 and 8.02 for Up = 1, 



Xmin 



Similarly, the 95.4% confidence level is determined by Ax = X 
2 and 3, respectively. In the MEF model, there are 2 free model parameters, namely ftmo and C- Note 
that E{z) for the MEF model has been given in Eq. ([251) . By minimizing the corresponding x^, we find 
the best-fit parameters Vlmo — 0.2704 and ( = 4 x 10""^, while Xmin — 542.683. In Fig.[Tl we present the 
corresponding 68.3% and 95.4% confidence level contours in the flmo — C parameter space. 

Next, we add the data from the observation of the large-scale structure (LSS). Here we use the distance 
parameter A of the measurement of the baryon acoustic oscillation (BAO) peak in the distribution of 
SDSS luminous red galaxies [2^, l30|, which contains the main information of the observations of LSS. 
The distance parameter A is given by 



A^nll^Eizt 



-1/3 



1 

Zb Jo 



dz 



E{z) 



2/3 



(34) 



8 



where zt = 0.35. In [ISl, the vahie of A has been determined to be 0.469 (ns/0.98)"°-3^± 0.017. Here the 
scalar spectral mdex Ug is taken to be 0.963, which has been updated from the WMAP 7-year (WMAP7) 
data [3i|. Now, the total x^ = X^ +xiss' '^here xfi is given in Eq. ([55]) . and Xlss ~ {■^~ -^obsY /cr\. By 
minimizing the corresponding x^, we find the best-fit parameters Vlmo — 0.2733 and C = 9 x 10^*, while 
Xmin — 542.734. In Fig. [2j we present the corresponding 68.3% and 95.4% confidence level contours in 
the J7mO — C parameter space. Comparing Fig. [5] with Fig. [1] it is easy to see that the constraints become 
much tighter. 

Then, we further add the data from the observation of the cosmic microwave background (CMB). 
Here we use the the shift parameter i?, which contains the main information of the observations of the 
CMB [alHll. The shift parameter R of the CMB is defined by [H,!!!] 



-ni/2 



R = Vl 



dz 

W) 



(35) 



where the redshift of recombination z* = 1091.3 which has been updated in the WMAP7 data |3l|. The 
shift parameter R relates the angular diameter distance to the last scattering surface, the comoving size 
of the sound horizon at z^. and the angular scale of the first acoustic peak in CMB power spectrum of 
temperature fluctuations fH,!!!]. The value of R has been updated to 1.725 ± 0.018 from the WMAP7 
data [m. Now, the total x"^ = xl + xiss + Xcmb^ where Xcmb = (-R ^ Robsf/o\. By minimizing the 
corresponding x^, we find the best-fit parameters Q.mo = 0.2699 and C, = 0.0165, while X^^^ — 542.879. In 
Fig.[3l we present the corresponding 68.3% and 95.4% confidence level contours in the ri„io — C, parameter 
space. Clearly, the constraints become tighter, and the best-fit C significantly deviates from zero. 



Vn 



0.2 


^ 


0.15 


/ \ 


O.I 


/ r'\ \ 


0.05 






' \\ 








• 





0.24 0.26 0.28 0.3 0.32 

■^mO 



FIG. 3: The same as in Fig. \T\ except that this result is obtained by using the combined observational data of 
557 Union2 SNIa, the distance parameter A from LSS, and the shift parameter R from CMB. 



Finally, we add the 59 Hymnium Gamma- Ray Bursts (GRBs) [28|, which can be used to constrain 
cosmological models without the circularity problem. In fact, GRBs are a complementary probe to SNIa 
(see e.g. [3J| and references therein), whose data points are also given in terms of the distance modulus 
l^obs{zi), similar to the case of SNIa. Therefore, the corresponding Xgrb from 59 Hymnium GRBs is 
also given in the same form of Eq. ([55]) . but the data points are replaced by the ones of GRBs. Now, 
the total X^ = xfi + Xlss + Xcmb + Xgrb- ^Y minimizing the corresponding x^, we find the best-fit 
parameters fimo — 0.2704 and ^ — 0.0188, while Xmin — 566.12. Note that the number of the total data 



points increased by 59, this Xmin is still very good. In Fig. |4l we present the corresponding 68.3% and 
95.4% confidence level contours in the ilmo — C parameter space. The difference between Figs. 2] and [3] is 
small, mainly due to the relatively weak constraint ability of current GRBs sample. This situation will 
be changed when more and more high-quality GRBs are available in the future. 
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FIG. 4: The same as in Fig. [TJ except that this result is obtained by using the combined observational data of 557 
Union2 SNIa, the distance parameter A from LSS, the shift parameter R from CMB, and 59 Hymnium GRBs. 



III. FURTHER DISCUSSIONS 



In the previous section, we have obtained the cosmological constraints on the MEF model. Here, we 
would like to continue with further discussions. 

One of the important issues is to justify the approximate solution E{z) given in Eq. (1281) . For the very 
small C, such as the C = 10~^ in |2y| or even smaller, we can of course safely use the approximate solution 



E{z) given in Eq. ([25]) . The question is when ( is not so small, can we still safely use the approximate 
solution E{z) given in Eq. ([^5])? As the first example, we consider C — 0.017 and flmo = 0.27, which is 
near to the best fits of the joint constraints from SNIa+LSS+CMB and SNIa+LSS+CMB+GRBs. We 
can numerically solve the exact differential equation (P5)) to find the exact solution of E{z) (of course, 
the initial condition is E{z = 0) = 1). To avoid confusion, we label the exact solution from Eq. (I?5|) and 
the approximate solution E{z) given in Eq. (^5)) as E^^ and Eapp, respectively. In Fig. [SJ we show the 
difference AE = Eapp — Eex and the relative difference AE/E = {Eapp — Eex)/Eex for the case with 
^ = 0.017 and ilmo = 0.27. Clearly, the difference between Eapp and Eg^ is very small, and hence we can 
reliably use the approximate solution E{z) given in Eq. (1^5)1 . Next, we consider a larger (. From Figs. [3] 
and[31 one can see that the upper edges of the 95.4% confidence level contours of the joint constraints 
from SNIa+LSS+CMB and SNIa+LSS+CMB+GRBs extend to C ~ 0.2. So, we choose the case with 
( = 0.2 and i^mo = 0.27 to be the second example. Again, we present the corresponding AE and AE/E 
in Fig. [SI One can see that the difference between Eapp and E^x is still fairly small, and hence we can 
also reliably use the approximate solution E{z) given in Eq. (|28p . Of course, if the computational ability 
of the computer becomes more powerful, it is undoubtedly the best choice to use the exact solution E^x 
from numerically solving Eq. (j23p . When we work with a less powerful computer, it is suitable to use the 
approximate solution E{z) given in Eq. (^5]) . 
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FIG. 5: The difference AE = Eapp — E^x and the relative difference AE/E = {Eapp — Eex)/Esx for the case with 
C = 0.017 and D,„ 



0.27. See the text for details. 



Next, we turn to the second issue. Although in the text below Eq. (j28p we have shown the key point 
to understand the reason for accelerating the universe without dark energy in the MEF model, it is still 
desirable to show the visualized plots. As is well known, the deceleration parameter is given by [3a| 



aiJ2 



H 



= -l-77^ = -l + (l 



z)E-' — 
dz 



(36) 



In Fig. [71 we plot the reduced Hubble parameter E{z) given in Eq. (125)) and the corresponding deceleration 
parameter q{z) for the cases with the best-fit parameters of the joint constraints from SNIa+LSS+CMB 
(black solid lines) and SNIa+LSS+CMB+GRBs (red dashed lines). In fact, the plot lines for both cases 
of SNIa+LSS+CMB and SNIa+LSS+CMB+GRBs are heavily overlapped. From Fig. [3 we can clearly 
see that the deceleration parameter q crosses the transition line g = at redshift Zt — 0.75, and the 
universe can really be accelerated in the late time without dark energy. By the way, we note that the 
usual relation 

WeS = -1 + 3(1 + ^)^"^ = 3 (^'^ " ^) 

does not hold in the MEF model, due to the fact that Friedmann equations have been modified. On the 
contrary, q{z) given in Eq. p6|) holds in any models since it comes from definition directly. 

Finally, as mentioned in Sec. HI the argument that MEF can avoid the conflict with the experiments 
testing the inverse square law relies on a small go- We know that go — Ho — C-^o- In |20|' ^^^ chose a 
tiny ^ = 10"^ for example, which can of course make a very small go- However, as shown in this work, 
the allowed C can be in the range < C^O-2 within 95.4% confidence level. Notice that the strength 
of gravitational fields are of order 10 N • kg~^ on the Earth, 10"^ N • kg~^ in the solar system |23- For 
the best-fit C ^ 10"^, the corresponding go = Ho = C-^o ^ 10"^^ N • kg~^. Even for the upper bound 
C ^ 10~^, the corresponding go = Ho = CHq ~ 10^^^ N • kg^^. Therefore, on the Earth or in the 
solar system, g 3> go, hence x = go/ g ^ 1, and then we have D[x) — >■ 1. Since the experiments testing 
the inverse square law were done on the Earth or in the solar system, the deviation from the inverse 
square law are extremely tiny. The significant deviation from the inverse square law can only occur in 
the very large scale in the universe where the gravitational fields are very weak, and hence it can escape 
the detection of the experiments testing the inverse square law which were done on the Earth or in the 
solar system. 
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FIG. 6: The same as in Fig. [S] except for the case with (" 
enlarged parts in the redshift range < z < 50. 
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IV. CONCLUSION AND REMARKS 



In summary, we considered the cosmological constraints on the MEF model in this work, by using the 
observational data of 557 Union2 SNIa, the distance parameter A from LSS, the shift parameter R from 
CMB, and 59 Hymnium GRBs. We found that the key parameter ( in MEF model has been limited in a 
narrow range < C^ 0-2 within 95.4% confidence level. By using the important result given in Eq. (05| . 
we have clearly shown the key point to understand the reason for accelerating the universe without dark 
energy in the MEF model. We showed the MEF model reduces to ACDM model when C <C 1. However, 
the best-fit C for the observations SNIa+LSS+CMB and SNIa+LSS+CMB+GRBs significantly deviates 
from zero. This indicates the new feature of MEF model different from ACDM model. We have justified 
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FIG. 7: The reduced Hubble parameter E{z) given in Eq. H28I1 and the corresponding deceleration parameter q{z) 
for the cases with the best-fit parameters of the joint constraints from SNIa+LSS+CMB (black solid lines) and 
SNIa+LSS+CMB+GRBs (red dashed lines). 



the approximate solution E{z) given in Eq. (P5| . We plotted E{z) and q{z) as functions of redshift z, 
and clearly showed that the universe can be accelerated in late time without dark energy. Finally, we 
have shown that MEF can avoid the conflict with the experiments testing the inverse square law. 

After all, some remarks are in order. In the MEF model, as shown in this work, the universe can 
be accelerated without dark energy. The only component is dust matter. The MEF model is in fact a 
modified gravity model, similar to the f{R) models and the braneworld models. The MEF model can be 
degenerate to ACDM model, but it has not an explicit cosmological constant in the model. This is a big 
advantage in fact, beyond some f{R) and braneworld models in this sense. 

Secondly, we point out the possibility to extend the original MEF model. In principle, it is not necessary 
to restrict the energy component in the universe to be dust matter only. The universe can contain other 
components, such as, dark energy. For example, we can consider a universe containing both dust matter 
and dark energy whose equation-of-state parameter (EoS) wx is a constant. In this case, the total energy 
density p = p,n + px ^ Pmo a~'^ + Pxo a~^^^^'^^\ Substituting into Eq. ((20|) . we have 



-2D 



3C/E 



'.C/E-1 



■2E 



dE 
dz 



= aanO (1 + Z)2 + 3 (1 + wx) f^X0(l + Z) 



3t«x+2 



(37) 



where 0,xo = {SttG pxo) / {^Hq) . Note that ftmo + ^xo 7^ Ij since Friedmann equation has been modified 
in the MEF model. In principle, one can numerically solve Eq. (1571) to get E{z). Similar to the original 
MEF model, we find that an approximation of the exact differential equation (I37p is given by 



1- -— ] ■2EdE ^n 

AE ' 



mO 



da 



-3 



nxoda-^'^^+'"''K 



Integrating Eq. 



we have 



E^--CE^ n,no a-^ + Qxo fl-^li+^-x) + const.. 



(38) 



(39) 



where const, is the integral constant, which can be determined by requiring E{z = 0) = 1. Finally, we 
find that 



E^-^CE^ n^o (1 + ^)' + ^xo (1 + z)3(i+--) + (i-^(- f]„o - n^^\ 



(40) 
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which is a quadratic equation of E in fact. Noting that E is positive, we solve Eq. (|40| to get 



E{z) = -C + 



1 (9 



I ^ C' + 4 n,no (1 + zf + fixo (1 + z)3(i+--) + (i-^(- r)„^o _ ri^o") | , 



2 4 



(41) 



In fact, this is just a simple example. One can include any type of dark energy, for instance, the CPL 
dark energy whose EoS is given by Wde = wa + Wa{l — a), quintessence, phantom, k-essence, hessence, 
(generalized) Chaplygin gas, holographic/agegraphic dark energy, vector-like dark energy, spinor dark 
energy, and so on. Therefore, we would like to give the more general formulae. In this case, the total 
energy density p = pm + Pde = Pmo a~'^ + Pde,o f{o), where /(a) can be any function of a which satisfies 
f{a— 1) = 1. Substituting into Eq. ^^, we obtain 



-2D 



;C/s~l 



dz 



nae,Q{l + z)-^f, 



(42) 



where /' = df /da, and ^defi = (87rGpde,o)/(3iJo). Note again that fimo + ^de,o 7^ Ij since Friedmann 
equation has been modified in the MEF model. Eq. (H^ is the exact differential equation, which can be 
used to find the exact E{z) numerically. Also, we give the corresponding approximate solution as 



^ ' 4^2 14^ 



f7,„o (1 + zf + Udefl I 



1 



1 



l< 



a 



rnO 



^de,0 



1/2 



(43) 



Similarly, if one need to add other component, such as radiation, it is not a hard work. In fact, we can 
give the most general formulae. In this case, the total energy density p = po /(a), where /(a) can be any 
function of a which satisfies /(a = 1) = 1. Substituting into Eq. (|20p . we have 



-2Dl^ 



3C/E 



iC/E-i 



2e'^ 
dz 



-noil 



'f: 



(44) 



where Q,q = (87rGpo)/(3-ffo)- Note again that J7o 7^ 1, since Friedmann equation has been modified in 
the MEF model, namely H^ ^ (87rGp)/3. Eq. (jH]) is the exact differential equation, which can be used 
to find the exact E{z) numerically. Also, we give the corresponding approximate solution as 



E{z) 



3^ 1 

4^+2 



C' + 4 



«o/ 



1 



1 



1 - ^ C - f^o 



1/2 



(45) 



In fact, noting that E — H/Hq, Eq. (HSI) can be regarded as the approximate modified Friedmann 
equation in the MEF model. If ^ <C 1, Eq. (HSI) reduces to 



^ = -IT- P + ^cS: 



(46) 



where Aeg — {1 — fto) Hq — const, is actually an effective cosmological constant. Therefore, in the 
most general case, we can clearly reveal the implicit root to accelerate the universe in the MEF model, 
regardless of the energy components in the universe. An effective cosmological constant is the intrinsic 
feature (we refer to e.g. [23| for a previous insight). The other exotic features of the MEF model could 
emerge only when (^ significantly deviates from zero. 

Thirdly, we would like to say some words on the understanding of the original entropic force model [l[ 
and the modified entropic force model ^2Q] . In fact, entropic force is just a new perspective to gravity, from 
the thermodynamical point of view. So, the original entropic force can only recover all results of the usual 
(Newton and Einstein) gravity. The only new thing is the reversed logic which might reveal the nature 
of gravity. In the original entropic force model [l|, using the fundamental assumptions Eqs. ((!]), ([3]), 
([5]) and ([6]) , Verlinde derived the Newton's law of gravitation Eq. ([8]) for the (non-relativistic) Euclidean 
spacetime in section 3 of [l|, and also derived the Einstein gravitational equations for any (relativistic) 
curved spacetime in section 5 of [l[ . On the other hand, the Friedmann equations were derived in [3, Q 
for the Friedmann-Robertson- Walker (FRW) universe. There is no any mixing here. We should mention 
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that both the original entropic force [l| and the modified entropic force [2^ cannot be understood in 
only non-relativistic or relativistic cases. In fact, they are equivalent to gravity itself in all cases. As the 
usual understanding, the Newton's law of gravitation is just the approximation of Einstein gravitational 
equations in the (non-relativistic) small scale limit, whereas the Friedmann equations are just the special 
case of Einstein gravitational equations in the cosmic scale (homogeneous and isotropic spacetime). The 
situation is similar in the modified entropic force model. Using the fundamental assumptions Eqs. ([T]), 
^ , ([5]) and pT|) , in 20] Gao derived the Newton's law of gravitation Eq. p^ for the (non-relativistic) 
Euclidean spacetime, and also derived the second Friedmann equation p^ for the FRW universe. Note 
that the first Friedmann equation can be derived from the second Friedmann equation (TTSl) and the 
energy conservation equation P^ . On the other hand, following Verlinde's derivations in section 5 
of [l| , one can derive the corresponding (modified) Einstein gravitational equations for any (relativistic) 
curved spacetime. In fact, this is just the lacked sector in the modified entropic force model. However, 
it is available in principle, although it has not been given in the literature. In the modified entropic 
force model, there is no any mixing too. The modified Newton's law of gravitation Eq. (jl4[) is just 
the approximation of the (lacked but available in principle) modified Einstein gravitational equations in 
the (non-relativistic) small scale limit, whereas the modified second Friedmann equation (J16p is just the 
special case of the (lacked but available in principle) modified Einstein gravitational equations in the 
cosmic scale (homogeneous and isotropic spacetime). 

Fourthly, we said that the MEF model is similar to /(i?)-gravity or braneworld scenario. Notice that 
they are similar only in the sense that the gravity has been modified in these models. Of course, both 
/(i?)-gravity and braneworld scenario were derived from the known actions, whereas the action for MEF 
is still lacked in the literature. However, as mentioned above, following Verlinde's derivations in section 5 
of [l|, in principle one can derive the corresponding (modified) Einstein gravitational equations for any 
(relativistic) curved spacetime. Once this lacked sector has been done, the explicit action is ready. Since 
the present work focuses on cosmology in the MEF model, we leave this task to future works. 

Fifthly, as mentioned in this work, in the MEF model, there is no dark energy in fact. The universe 
is matter-dominated always. The expansion of our universe is accelerated due to the fact that gravity 
has been modified. In the non-relativistic case, there is no dark energy too, but gravity is also modified. 
However, as mentioned in this work, this modification to Newtonian gravity is negligible on the Earth or 



in the solar system. In the larger scale, the modified gravity is described by Eq. (|14p . As shown in [37 1 
the Debye entropic force can be an alternative to the modified Newtonian dynamics (MOND) to explain 
the rotational velocity curves of spiral galaxies. In fact, the MEF model '2(| and the Debye entropic force 
model [33] are very similar. So, it is anticipated that the "non-relativistic cosmology" of the MEF model 
could be an alternative to dark matter, which is usually invoked to explain the rotational velocity curves 
of spiral galaxies. 

Finally, we admit that the entropic force proposed by Verlinde is based on several unproved hypotheses, 
and it is still controversial in the physical community. On the other hand, the Debye model in the 
thermodynamics has not been used in the gravity theory previously. However, in the history, many great 
theories also appeared controversially in their beginning. Therefore, we consider that it is better to keep 
an open mind to these speculative attempts. 
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